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Abstract 

Consider a bounded planar domain D, an instance h of the Gaussian free field 
on D, with Dirichlet energy (2ir)^ 1 f D Vh(z) ■ Vh(z)dz, and a constant < 7 < 2. 
The Liouville quantum gravity measure on D is the weak limit as e —* of the 

measures 

where dz is Lebesgue measure on D and h E (z) denotes the mean value of h on the circle 
of radius e centered at z. Given a random (or deterministic) subset X of D one can de- 
fine the scaling dimension of X using either Lebesgue measure or this random measure. 
We derive a general quadratic relation between these two dimensions, which we view 
as a probabilistic formulation of the KPZ relation from conformal field theory. We also 
present a boundary analog of KPZ (for subsets of dD). We illustrate (via heuristics and 
announced results) the connection between discrete and continuum quantum gravity 
and provide a framework for understanding Euclidean scaling exponents via quantum 
gravity. 
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"There are methods and formulae in science, which serve as master-keys to 
many apparently different problems. The resources of such things have to be 
refilled from time to time. In my opinion at the present time we have to develop an 
art of handling sums over random surfaces. These sums replace the old-fashioned 
(and extremely useful) sums over random paths. The replacement is necessary, 
because today gauge invariance plays the central role in physics. Elementary 
excitations in gauge theories are formed by the flux lines (closed in the absence 
of charges) and the time development of these lines forms the world surfaces. All 
transition amplitude[s] are given by the sums over all possible surfaces with fixed 
boundary." (A.M. Polyakov, Moscow, 1981.) [PolSlaj 



1 Introduction 



1.1 Overview 



The study of certain natural probability measures on the space of two dimensional Rie- 
mannian manifolds (and singular limits of these manifolds) is often called "two-dimensional 
quantum gravity." These models have been very thoroughly studied in the physics liter- 
ature, in part because of connections to string theory and conformal field theory [Pol81a 



IPolSlbl lPol87al lPo!89l ISeTQOl IGM931 IDa^94l IDav951 IA.TW951 IAW951 IDFGZ.T951 lKle95l IKH961 
IADJ97[ EynOl , Dup 06j , and to random matrix theory and geometrical models; see, e.g., the 



references |BIPZ78llADF85l IKKM851 IDav85l IBKKM86al IBKKM86bl IKaz86l IDK88al IDK901 
IGK891 IKos89a| IKos89bj IMSS91L lKK92j IEZ921 IDau951 lEK95j lKH95j IBDKS951 IAAMT96} 



Dup98| |Dup99a| |Dup99b[ |Dup99c[ [EB991 [KZJ99J EqsQQI |DupOO| IDFGGOO} |Dup04| . More 



recently, a purely combinatorial approach to discretized quantum gravity has been success- 
ful |Sch98l IBS021 IBDFG021 IA"S03l IBDFG03al IBDFG03bl IDFG051 IBDFG071 lMie07l ILG071 
IMM071 lBG08al IMW081 lMie08l IBG08bl ILG08] . 

One of the most influential papers in this field is a 1988 work of Knizhnik, Polyakov, and 
Zamolodchikov [KPZ88J. Building on a 1987 work of Polyakov [Pol87b], the authors derive a 
relationship (the KPZ formula) between scaling dimensions of fields defined using Euclidean 
geometry and analogous dimensions defined via Liouville quantum gravity (as described 
earlier in [Pol8 1 at IP0I8 1 b j ) . An alternative heuristic derivation using Liouville field theory in 
the so-called conformal gauge was proposed shortly after |Dav88at IDK89] (see also |Tak93j ). 
The original work by KPZ has been cited roughly a thousand times in a variety of contexts, 
which we will not attempt to survey here, though we mention that there have been a number 
of explicit calculations in Liouville field theory with matching results in the random matrix 
theory approach, e.g., [DO94llZ^II^lTisTmiPT02llK^^ 

The relationship in [KPZ 88] has never been proved or even precisely formulated mathe- 
matically. The main goal of this work is to formulate and prove the KPZ scaling dimension 
relationship in a probabilistic setting. 

1.2 Critical Liouville quantum gravity 

The study of two dimensional random surfaces makes frequent use of the Riemann uni- 
formization theorem, which states that every smooth simply connected Riemannian mani- 
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fold Ai can be conformally mapped to either the unit disc D, the complex plane C, or the 
complex sphere C U {oo}. (If a manifold is not simply connected then its universal cover 
can be conformally mapped to one of these spaces. See, e.g., Chapter 4 of |FK92] for more 
exposition; see also [WGY051 IJWGY051 IGWY031 IGY021 IDL,I+07| for approximation algo- 
rithms and beautiful computer illustrations of these maps.) Another way to say this is that 
Ai can be parameterized by points z = x + iy in one of these spaces in such a way that the 
metric takes the form e x ^(dx 2 + dy 2 ) for some real- valued function A. The (x, y) are called 
isothermal coordinates or isothermal parameters for Ai. In most of this paper we let the 
parameter space be a general simply connected proper subdomain D of the plane (which, of 
course, is conformally equivalent to D). 

We remark that the existence of isothermal coordinates does not require that Ai be 
smooth; for example, it can be deduced whenever Ai can be parameterized by a simply 
connected planar domain in which the metric has the form E(x,y)dx 2 + 2F(x,y)dxdy + 
G(x,y)dy 2 where EG — F 2 > 0, E > 0, and E, F, and G are /3-Holder continuous for some 
< /3 < 1 |Che55j . 

Length, area, and curvature are easy to express in isothermal coordinates. The length of 
a path in Ai parameterized by a smooth path P in D is given by 

J e x ^ 2 ds, 

where ds is the Euclidean length measure on D. Given a measurable subset A of D, the 
integral J A e x ^dz (where dz denotes Lebesgue measure on D) is the area of the portion of 
Ai parameterized by A. The function K = — e~ A AA (where AA = \ xx + \ yy is the Laplacian 
operator) is called the Gaussian curvature of Ai. If A is a measurable subset of the (x, y) 
parameter space, then the integral of the Gaussian curvature with respect to the portion 
of Ai parameterized by A can be written J A e x ^K{z)dz = J A —AX(z)dz where dz denotes 
Lebesgue measure on D. In other words, — AA gives the density of Gaussian curvature in 
the isothermal coordinate space. In particular, Ai is flat if and only if A is harmonic. 

The above suggests that one can study random simply connected Riemannian manifolds 
by studying random functions A on C or C U {oo} or any fixed simply connected subdomain 
D of C. In the probabilistic formulation of the so-called critical Liouville quantum gravity, A 
is taken to be a multiple of the Gaussian free field (GFF), although some care will be required 
to make sense of this construction, since the GFF is a distribution and not a function. (The 
relationship between our probabilistic formulation and the original formulation of Polyakov 
will be discussed in Section EJ) 

For concreteness, let h be an instance of a centered GFF on a bounded simply connected 
domain D with zero boundary conditions. This means that h = Yl n a nfn where the a n are 
i.i.d. zero mean unit variance normal random variables and the f n are an orthonormal basis, 
with respect to the inner product 

(/ 1 ,/ 2 )v:=(27r)- 1 / Vf 1 (z)-Vf 2 (z)dz, 

J D 

of the Hilbert space closure H(D) of the space H S (D) of C°° real- valued functions compactly 
supported on D. Although this sum diverges pointwise almost surely, it does converge almost 
surely in the space of distributions on D, and one can also make sense of the mean value 
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of h on various sets. (See |She07] for a detailed account of this construction of the GFF; 
see Section [3TT1 for a quick overview. Note that the (27r) _1 in the definition above does not 
appear, e.g., in [Shc07j; including this factor in the definition, as is common in the physics 
literature, is equivalent to multiplying the corresponding h by ^/2n. This will simplify some 
of our formulas later on. In particular, in this formulation the two point covariance scales 
like — log(|z — w\) instead of — (27r) _1 log(|z — w\); see Section [37TI ) 

Given an instance h of the Gaussian free field on D, let h e (z) denote the mean value of h 
on the circle of radius e centered at z (where h(z) is defined to be zero for z £ C\D). This is 
almost surely a locally Holder continuous function of (e, z) on (0, oo) x D (see Section I3~TT) . 
For each fixed e, consider the surface Ai e parameterized by D with metric e^^' (dx 2 + dy 2 ) . 
We would like to define a surface Ai parameterized by D to be some sort of limit as e — > 
of these surfaces. Since we would not expect the limit to be a Riemannian manifold in any 
classical sense, we have to state carefully what we mean by this. There are many ways we 
could attempt to make sense of this limit, depending on what quantities we focus on. For 
example, we could consider 

1. The length of the shortest path connecting a fixed pair of points in D. 

2. The area of a fixed subset of D. 

3. The length of a fixed smooth curve in D. 

4. The length of a smooth boundary arc of D (which becomes interesting when h is an 
instance of the GFF with free boundary conditions). 

Intuitively, we might expect each quantity above to scale like a random constant times a 
(possibly different) power of e as e tends to zero — i.e., we would expect that if the A4 e 
were rescaled by the appropriate powers of e, the above quantities would have limits as 
e — > 0. Focusing on lengths of shortest paths, one might guess that the random surfaces A4 £ 
(rescaled by some power of e) would almost surely converge (in some natural topology on the 
set of metric spaces) to a non-trivial random metric space parameterized by D. However, 
this is not something we are currently able to prove. Focusing on areas, one might expect 
that for some a the renormalized area measures e a e lhe ^dz would almost surely converge 
weakly to a random measure on D. This is the limit we will construct and work with in this 
paper. We will also address the lengths of fixed curves and boundary curves; see Section 
[6j Although the constructions are quite similar, we will not use the so-called Wick normal 
ordering terminology in this paper (see e.g., [Sim74] ) . We present a self-contained proof 
of the following (although similar measures have appeared much earlier, and are called the 
H0egh-Krohn model |HK71] - see also [AGHK79} IAHK74] for a discussion on the level of 
Schwinger functions, and a more recent survey [AHKPS92J): 

Proposition 1.1. Fix 7 £ [0,2) and define h and D as above. Then it is almost surely the 
case that as e — > along powers of two, the measures jj £ := e 7 l 2 e~ <he ^dz converge weakly to 
a limiting measure, which we denote by jj = fih = e lh ^dz. This remains true if we replace 
h with a non-centered GFF on D — i.e., if we set h = h + h° where h is the zero boundary 
GFF on D and h° is a deterministic, non-zero continuous function on D. 
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For each z G D, denote by C(z; D) the conformal radius of D viewed from z. That is, 
C(z;D) = where <fi : D — > D is a conformal map to the unit disc with </>(z) = 0. 

The following gives an equivalent definition of \i. 

Proposition 1.2. Write h = h + h° where h is the zero boundary GFF on D and h° is a 

deterministic continuous function on D. Let /i,/2, . . . be an orthonormal basis for H(D) 
comprised of continuous functions on D and let h n be the expectation of h given its projection 
onto the span of {fx, / 2 , . . . , f n }. (In other words, h n is h° plus the projection of h onto the 
span of {fi, f 2 , ■ ■ ■ , f n }-) Then \i = Hh (as defined in Proposition !! . 1\) is almost surely the 
weak limit for n — > +00 of the measures 

fx n = exp (^h n (z) - -^V&rh n (z) + ^ logC(z; D^J dz. (1) 

For each measurable A C D , we have K[fi(A)\h n ] = fi n (A). In particular, 

Efx(A) = [ C(z;D)^e^ h " (z) dz. 

J A 

Intuitively, we interpret the pair (D,fi) as describing a "random surface" Ai parameter- 
ized conformally by D, with area measure given by /i. The term "random metric" is often 
used as well; however, we stress that, since we have not endowed D with a two point distance 
function, "random metric" in the Liouville quantum gravity context does not mean "random 
metric space." 

1.3 Scaling exponents and KPZ 

Definition 1.3. For any fixed measure fi on D (which we call the "quantum" measure), we 
let B s (z) be the Euclidean ball centered at z whose radius is chosen so that fj,(B s (z)) = 5. (If 
there does not exist a unique 5 with this property, take the radius to be sup{e : fi(B e (z)) < 5}.) 
We refer to B s (z) as the isothermal quantum ball of area 6 centered at z. In particular, 
if j = then fi is Lebesgue measure and B 5 (z) is B e (z) where 5 = ne 2 . 

Given a subset X C D, we denote the e neighborhood of X by 

B £ (X) = {z:B £ (z)nX^®}. 

We also define the isothermal quantum 5 neighborhood of X by 

B 5 (X) = {z : B 5 (z)nX^ 0}. 

Translated into probability language, the so-called KPZ formula is a quadratic relation- 
ship between the expectation fractal dimension of a random subset of D defined in terms of 
Euclidean measure (which is the Liouville gravity measure with 7 = 0) and the corresponding 
expectation fractal dimension of X defined in terms of Liouville gravity with 7 7^ 0. 

Fix 7 G [0, 2) and let denote Lebesgue measure on D. We say that a (deterministic 
or random) fractal subset X of D has Euclidean expectation dimension 2 — 2x and 
Euclidean scaling exponent x if the expected area of B e (X) decays like e 2x = (e 2 ) x , i.e., 

e^o lose 2 
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We say that X has quantum scaling exponent A if when X and // (as defined above) are 
chosen independently we have 

Um logEMS'W) = A 
<5^o log 5 

(Section [7] will provide some discrete quantum gravity heuristics that motivate the idea of 
taking X and D independent of one another, as well as our particular definition of scaling 
exponent.) 

The following is the KPZ scaling exponent relation. To avoid boundary technicalities, 
we restrict attention here to a compact subset of D. The case of boundary exponents will 
be dealt with in Section O 

Theorem 1.4. Fix 7 G [0, 2) and a compact subset D of D. If X n D has Euclidean scaling 
exponent x > then it has quantum scaling exponent A, where A is the non-negative solution 
to 

7\ 2 , f-, 7 2 



* = T A^l- T jA. (2) 

It also turns out that Theorem 11.41 admits the following straightforward generalization: 

Theorem 1.5. Let X be any random measurable subset of the set of all balls of the form 
B E (z) for e > and z in a fixed compact subset D of D. Fix 7 G [0, 2). Then if 

e^o loge z 

then it follows that, when X and fj, (as defined above) are chosen independently, we have 

,. logEyujz : B 5 (z) G X} 

hm Pl = A, 

s^o log 5 

where A is the non-negative solution to 



.-^■ + (l-f)A. 

(Note that expectation in the above theorem is with respect to both random variables, 
X and /x.) We obtain Theorem 11.41 as a special case of Theorem 11.51 by writing X = {B £ (z) : 
B £ (z) nl^ 0}. Theorem 11.51 allows us to consider x that are greater than 1 (in which case 
the "dimension" 2 — 2x would be negative). If one considers, for example, a conformal loop 
ensemble on D with k = 6 (corresponding to a scaling limit of the cluster-boundary loops 
in site percolation on the triangular lattice) one could let X be the set of balls contained 
in D that intersect i distinct "macroscopic" loops (where "macroscopic" means that their 
diameters are greater than some fixed constant). In this case, the value x depends on I 
and is called a multi-arm exponent [SD87, Dup99a, A DA991 ISWOlj and we may view the 



corresponding A as a quantum analog of such an exponent. 

As another example, for some integer L fix distinct points zi, Z2, ■ ■ ■ , Zi in D\ D and run 
L independent Brownian motions started at the points Z\, . . . , Zl. Then let X be the set of 
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balls B E (z) contained in D with the property that the Brownian motions — stopped at the 
first time they intersect dB E (z) — do not intersect one another. 

In this case, the Euclidean scaling exponent x = xl is called a Brownian intersection 
exponent intersection. It was conjectured in [DK88b| and rigorously derived in a celebrated 
series of papers by Lawler, Schramm, and Werner using the Schramm-Loewner evolution with 
k = 6 jLSWOlal ILSWOlbl ILSW02] : 

*l = ^(4L 2 -1). 

Although we will not fully explain this in this paper, there is a close connection between 
SLE re and Liouville quantum gravity models with 7 = ->/min{fi;, 16/k} (see Section [8]), in 
agreement with the relationship between CFT central charge c and parameter 7 in Liouville 
quantum gravity [KPZ881 IDav88bl IDK891 1551501 IGM93j . Taking 7 = = ^/8/3 and 

as above, the KPZ formula gives 

2 V 2 

which is an affine function of L. The first co-author predicted several years ago, based on an 
approach via discrete quantum gravity models, that this A would be an affine function of L 
(see |Dup98 , Dup99b Dup99c Dup04| and the discussion in Section [7]). The derivation is 



based on a simple and general geometric argument that discrete quantum gravity exponents 
should be in a certain sense additive together with a heuristic connection between the discrete 
and the continuous models. A direct calculation via discrete graphs appears in |Dup98|. 
This is related to the cascade relations given earlier by Lawler and Werner using different 
techniques [LW99j . 

Two papers that build on our work (as announced and presented in talks and minicourses 
beginning in 2007) have already been posted online: Benjamini and Schramm cited the ideas 
of our paper to produce an analog of Theorem 11.41 in a one dimensional cascade model; their 
proof uses a Frostman measure construction in place of the large deviations construction 
used here, and almost sure Hausdorff dimension in place of expectation dimension [BS08] . A 
follow up paper [RV08] adapts the arguments of [BS08] to a class of cascade models, which 
was expanded to include (in a revised version) a measure based on the exponential of the 
Gaussian free field, like the measures we construct here. 

Intuitively, one reason to expect Hausdorff-like variants of KPZ to be accessible is that 
the second moments (and higher moments) of the random measures are essentially trivial to 
compute (see Section 1372]) . It might be interesting to try to derive other variants of KPZ: 
for example, one could try to relate the actual Minksowki or Hausdorff measure of a set, 
in the Euclidean sense, with some kind of expected Minkowski or Hausdorff measure in the 
quantum sense. We will not address these alternative formulations here. However, we will 
present below a picturesque formulation of KPZ in terms of box decompositions. 

1.4 Statement of box formulation of KPZ 

Define a diadic square to be a closed square (including its interior) of one of the grids 
2 _fc Z 2 for some integer k. Let /1 be any measure on C. For 5 > 0, we define a (fi, 5) box 
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S to be a diadic square S with /jl(S) < 5 and /i(S") > 5 where S' is the diadic parent of 
S. Clearly, if a point z G C does not lie on a boundary of a diadic square — and it satisfies 
fJ-({z}) < 5 < /i(C) — then there is a unique (/i, 5) box containing z, which we denote by 
S\z). Let be the set of all (fi,5) boxes. The boxes in do not overlap one another 
except at their boundaries. Thus, they form a tiling of M 2 (see Figures [TJ, [2} and [3] for an 
illustration of this construction on a torus). 

We remark that the (/i, 5) boxes should not be confused with the diadic boxes in the so- 
called 5-Calderon Zygmund decomposition of \i. Readers familiar with that decomposition 
may recall that while the (//, 5) boxes are diadic squares S with fi(S) < 5 < /x(S"), the 8- 
Calderon Zygmund boxes are diadic squares S with fj,(S) / f/,o(S) > 5 > fi(S') / ^(S 1 ) , where 
Ho is Lebesgue measure. Roughly speaking, the \i measure on each ([/,, S) box approximates 
S, while the \x density on each Calderon Zygmund box approximates 5. 

When e is a power of 2, analogously define S £ (z) to be the diadic square containing z 
with edge length e. Likewise, define 

S £ (X) = {z:S £ (z)nX^®}, 

S 5 (X) = {z s s (z)nx ^®}. 

The following gives the equivalence of the scaling dimension definition when boxes are used 
instead of balls. (The first half is well known and easy to verify.) 

Proposition 1.6. Fix 7 G [0,2) and let X be a random subset of a deterministic compact 
subset of D. Let N(fi,5,X) be the number of (n,5) boxes intersected by X and N(e,X) the 
number of diadic squares intersecting X that have edge length e (a power of 2). Then X has 
Euclidean scaling exponent x > if and only if 

^logEp^X)] 
e^o log e 2 



or equivalently, 



ioge 2 



£ 



Similarly, X has quantum scaling exponent A if and only if when X and fi (as 
defined above) are chosen independently we have 



or equivalently, 



]im ^nsN( flS ,x)] = 

6-+0 log 5 

lk iD 8 E[%,q)] _ A 1 

-5^0 log 5 

Of course, this immediately implies the following restatement of Theorem 11.41 in terms of 
boxes instead of balls: 

Corollary 1.7. Fix 7 G [0, 2) and a compact subset D of D and X and \x as above. Then if 

, im i°ggvM01 = x _ L 



e- 



+0 loge 2 
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Figure 1: (//, 5) boxes of the random measure \x = e lh dz, where 7 = .5 and h is the (discrete) 
Gaussian free field on a very fine (1024 x 1024) grid on the torus, dz is counting measure 
on the vertices of that grid, and 5 is 2~ 12 times the total mass of /1. (We view \i as an 
approximation of the continuum Liouville quantum gravity measure.) One way to construct 
this figure is to view the entire torus as a square; then subdivide each square whose \i measure 
is at least 5 into four smaller squares, and repeat until all squares have \i measure less than 
5. The squares shown have roughly the same /1 size — in the sense that each square has fi 
measure less than 5 but each square's diadic parent has \i measure greater than 5. 
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Figure 2: Analog of Figure [U with 7 = 1, using the same instance h of the GFF. 
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Figure 3: Analog of Figure CD with 7 = 1.5, using the same instance h of the GFF. 



11 



for some x > then 



<5^0 log 5 

where A is the non-negative solution to (TJ|). 

One could also phrase Theorem 11.51 in terms of boxes instead of balls, but for simplicity 
we will refrain from doing this here. 



2 Coordinate changes and the physical Liouville action 

Polyakov understood early on that the Liouville quantum gravity action becomes a free field 
action in the conformal gauge, but he did not construct the random metric the way we do. 
In |Pol87bj . where Polyakov begins the KPZ derivation, he refers to the Liouville quantum 
gravity action and writes 

"The most simple form this formula takes is in the conformal gauge, where 
9ab = e^5 a fc where it becomes a free field action. Unfortunately this simplicity is 
an illusion. We have to set a cut-off in quantizing this theory, such that it is 
compatible with general covariance. Generally, it is not clear how to do this. For 
that reason, we take a different approach. " 

Indeed, the actual derivation given in [Pol87bJ and subsequently in Knizhnik, Polyakov, and 
Zamolodchikov [KPZ88J is more complicated than ours and is not based on the Gaussian 
free field. It does not give precise mathematical meaning to the random surfaces. We feel 
that the Gaussian free field based random metric we construct is the correct one, at least in 
the sense that it is likely to arise as a scaling limit of the discrete quantum gravity models 
mentioned in |KPZ88j (see Section [7j). In a way our approach is more similar to the work 
of David |Dav88aj and of Distler and Kawai [DK89], which heuristically derived KPZ from 
Liouville field theory in the so-called conformal gauge. 

In this section, we describe how the Liouville quantum gravity measure we construct 
transforms covariantly under coordinate changes and use this to explain the connection 
between the Gaussian free field and the more familiar and more general curvature-based 
definition of the Liouville action that is conventional in the physics literature. The covariance 
properties of the random metrics in our point of view are very simple and agree with those 
postulated in the physics literature. 

If is a conformal map from D to a domain D and h is a distribution on D, then we 
define the pullback h o <p~ l of h to be a distribution on D defined by (h o ml ,p) = (h, p) 
whenever p e H S (D) and p = \<f>'\~ 2 p o (Here <f> is the complex derivative of 0, and 
(h, p) is the value of the distribution h integrated against p.) Note that if h is a continuous 
function (viewed as a distribution via the map p — ► f D p(z)h(z)dz), then the distribution 
h o _1 thus defined is the ordinary composition of h and (viewed as a distribution). 

The following transformation rule is a simple consequence of Proposition 11.21 and the 
definitions above. 



12 



Proposition 2.1. Let h be an instance of the GFF on D and ip a conformal map from a 
domain D to D . Write h for the distribution on D given by h o ip + Qlog \ip'\ where 

2 7 

7 2 

Then fih is almost surely the image under ip of the measure \i~ h on D. That is, tt^(A) = 
Hh{i/j{A)) for each Borel measurable A C D. 



Proof. Using the notation of Proposition [L2l if /i, /2, • • • are an orthonormal basis for H(D), 
then the conformal invariance of (■, -)y implies that f\oip, f2°ip, ■ ■ ■ are an orthonormal basis 
for H(D), and as n — > oo the functions h n o ip converge in law to the GFF on D, and the 
functions h n = h n o ip + Qlog^'l converge in law to h. If we define ft n analogously to fi n 
in ([1]) but with h n replaced by h n , then the jx n converge weakly to the random distribution 

To see that \i is the image of \x under ip, we will observe that fi n is the image of [i n 
under ip for each n. To see this, consider the term Q log \ip'\ = (2/7) log \ip'\ + (7/2) log \ip'\ 
in the definition of h. Adding (2/7) log \ip' \ to h n o ip corresponds to multiplying (pQ) by 
a factor of \ip'\ 2 . This compensates for the fact that the Radon-Nikodym derivative of 
a measure on D at a point z and the Radon-Nikodym derivative of the same measure 
pushed forward on D at ip(z) differ by a factor of |^'(2;)| 2 . Adding ( / ~f / 2) log \ip'(z)\ to 
h n o ip compensates the expression (CEJ) for the change in conformal radius: log C(ip(z); D) — 
\ogC(z-D) = \og\iP\z)\. □ 

We interpret Proposition 12.11 as a rule for changing the parametrization of a random 
metric. For example, consider the random metric one gets by taking D to be some fixed 
domain. Then if we are given any other domain D and a conformal map ip : D — > D, we 
may wish to consider the same random metric parameterized by D instead of D. In this 
case, the transformation rule tells us that on D we should consider the Liouville quantum 
gravity measure defined using h = hoip + Q log \ip'\, where h o ip is the GFF on D with zero 
boundary conditions. 

We remark that one can make a similar argument when D is a curved simply connected 
manifold and ip : D — > D a conformal map; in this case, the function A can be defined 
as follows: if the metric 011 D — when parameterized by D using the map ip^ 1 — takes 
the form e xi - z \dx 2 + dy 2 ), for z E D, then we write X(ip(w)) = —2\og\ip'(w)\ for w E D. 
Although we will not prove it here, the analog of Proposition II . II for smooth curved surfaces 
is straightforward, and the transformation rule Proposition 12.11 remains the same in this case; 
as in the flat case, the law of the Liouville quantum gravity measure on D pulled back to D 
is that of h = hoip + Q log \ ip'\ where ho ip is the GFF on D with zero boundary conditions. 
(Alternatively, we may take this as a definition of the Liouville quantum gravity measure on 
curved D with zero boundary conditions.) 

The remainder of this subsection describes the connection between our notation and the 
common physics literature Liouville gravity notation. (This discussion can be skipped, on a 
first read, by readers with no prior familiarity with the latter.) What we call the GFF on D 
(with the 1/27T normalization, as discussed in the introduction) is often written (sometimes 
without a rigorous definition) as the measure e~ s ^dh, where 

S(h) = J Vh{z) ■ Vh{z)dz 
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is called the action and dh is denned heuristically as a "uniform measure on the space 
of all functions." (Of course, the latter makes perfect sense if one considers only a finite 
dimensional vector space of functions, such as real-valued functions defined on the vertices 
of a lattice, or functions whose Fourier coefficients beyond a certain frequency threshold are 
identically zero — in this case dh would be the Lebesgue measure on the vector space.) In 
this paper, we will write 

(h, h) v := — [ Vh(z) ■ Vh(z)dz, 
2vr J D 

so that the above becomes S(h) = \{h, h)\j. 

In the following, let D be a subdomain of C and D a possibly curved surface for which 
there is a conformal map ip : D — > D. Write h° = log Now, if we switch parametrization 
to D, we are adding Qh° deterministically to h o to get h, so we may rewrite the action as 

S= ^(h-Qh°,h-Q~h°) v , 

which (at least when h° is smooth and compactly supported) is seen by integrating by parts 
to be equivalent (up to the additive constant |||<5^°||v) to 

S = -L Jdw (Vh(w) ■ Vh(w) + 2h(w)QAh°(wfj , (4) 

where the pairing Vh(w)-Vh(w) and the Laplacian Ah°(w) are now defined using the metric 
on D and where now dw represents the measure on D instead of D. This can also be written 

S = l- J dw (yh(w) ■ Vh(w) + Qh{w)K{w)^ , (5) 

where K is the Gaussian curvature of D and dw is integration with respect to that metric. 
(When h° is not compactly supported, the formula can be modified to include a term for 
boundary curvature, but we will not discuss this here.) 

Adding in one additional term which is a constant \x times the total area of D (and 
making the following symbol substitutions: b = 7/2, = h, g is the underlying metric of 
D, and j and k are summed-over indices ranging over the two tangent space directions), we 
obtain the more familiar formula for the Liouville action: 

S = ~L In ^ jk d^d k ip + QK V + 47r/xe 2 ^) , 

where Q = b + b~ l . The action is defined similarly when free boundary conditions are used 
instead of zero boundary conditions — or when D is a compact Riemann surface of some 
genus. (In this case, e" is an infinite measure, although it can be "localized," e.g., by 

requiring the mean value of tp to be zero.) 

This paper will focus exclusively on the case 7 G [0, 2) (which is said to correspond to 
physical models below the central charge c = 1 threshold) and fi = (the so-called critical 
Liouville quantum gravity). The string theory and quantum gravity literatures deal with 
other parameter choices as well — including non-zero ji and complex values for 7 and Q 
- but these appear to be beyond the scope of our methodology, in part because, when 
S is complex valued, the expression e~ s ^dLp is no longer a probability measure in even a 
heuristic sense. 
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3 Constructing the random measures 



3.1 GFF definition and normalization 

Let D be a bounded planar domain and let dz denote Lebesgue measure on D. We assume 
the reader is familiar with the Gaussian free field, as defined, e.g., in |She07] . but we briefly 
review the definition here. As described earlier, to make our formulas consistent with the 
physics literature, the definitions of Green's function and the Dirichlet form will differ from 
the ones in |She07j by factors of 2tt. 

Let H S (D) be the space of C°° real-valued functions compactly supported on D. We 
define the Dirichlet inner product 

(/i,/a)v := (27T)- 1 f Vft(z)-Vf 2 (z)dz, 

J D 

on H S (D). Then an instance h of the Gaussian free field (GFF) may be viewed as a standard 
Gaussian on the Hilbert space closure H(D) of H S (D) (i.e., as a sum of the form ^ n a n f n 
where /„ are any orthonormal basis for H(D)) — the sum converges almost surely in the 
space of distributions on D, see |She07] . In fact, we may define (h, f)y as random variables 
for non-smooth / as well; these are zero mean Gaussian random variables for each / G H(D), 
and 

Cov((/i,A)v,(A,/ 2 ) v ) = (/i,/ 2 )v. 

The collection of random variables (h, /)y for / G H(D) is thus a Hilbert space (isomorphic 
to H(D)) under the covariance inner product. 

When x G D is fixed, we let G x (y) be the harmonic extension to D of the function on 
dD given by — log \y — x\. Then Green's function in the domain D is defined by 

G(x, y) = - log \y-x\- G x (y). 

When x G D is fixed, Green's function may be viewed as a distributional solution of 
AG(x,-) = —2n5 x (-) with zero boundary conditions |She07] . It is non-negative for all 
x,y G D. 

For any function p on H S (D), we write A _1 p for the function 

-7T- / G{-,y) p{y)dy. 
^ Jd 

Integration by parts shows that this is a C°° (though not necessarily compactly supported) 
function in D whose Laplacian is p. We use the same notation for more general measurable 
functions p, as well as the case that p is a measure. (For example, we will sometimes speak 
of the inverse Laplacian of uniform measure on a particular circle or disc contained in D.) 

If fi = — A _1 pi and fi = — A _1 p 2 , then integration by parts implies that (/i,/2)v = 
(27r) _1 (pi, — A -1 /^), where (•, •) denotes the standard inner product on L 2 (D). We next 
observe that every h G H(D) is naturally a distribution, since we may define the map (h, •) 
by (h, p) := 2n(h, — A~ 1 p)y. (It is not hard to see that — A _1 p G H(D), since its Dirichlet 
energy is given explicitly by (jBD.) When —Af = p, we may write (h,p) = 2n(h,f)^, and 
hence 

Cov((h, Pl ),(h,p 2 )) =(27r) 2 (A,/ 2 ) v . 
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We claim that the latter expression may be rewritten to give 

Cov((/i,pi), (h,p 2 )) = / p 1 (x)G(x,y)p 2 (y) dxdy (6) 

JDxD 

where G(x,y) is Green's function in D. Since AG(x, •) = — 2tt5 x (-) and 

/ G(x,y) p 2 (y)dy = -2irA' 1 p 2 (x), 

J D 

we obtain ([6]) by multiplying each side by — Afi(x) = pi(x) and integrating by parts with 
respect to x. 

Denote by h £ (z) the average value of h on the circle of radius e centered at z. Similar 
averages were considered in |Bau90j . (For this definition, we assume h is identically zero 
outside of D.) Then h £ (z) is a Gaussian process with covariances defined by 

G £1 , £2 (zi, z 2 ) := Cov^Ox),/^^)) (7) 

given by 

J G{x,y)p z e \{x)p z £ l{y)dxdy 

where p z (x)dx is the uniform measure (of total mass one) on dB £ (z). In fact (like Brownian 
motion) the process h £ (z) determines a random continuous function (of z and e): 

Proposition 3.1. The process h £ (z) has a modification which is almost surely locally r\- 
Holder continuous in the pair (z,e) G C x (0, oo) for every rj < 1/2. 

In other words, the Holder regularity enjoyed by h £ ( function of the pair (e, z) - 

is the same as that of Brownian motion or the Brownian sheet. In fact, as we observe below 
(Proposition 13. 3p . when z is fixed, h e is a Brownian motion with respect to the parameter 
t = — logs. We may view h e as an approximation to h that gets better as e — > 0. Before 
we prove Proposition 13.11 let us make some observations about the covariance function 
G £lt£2 (zi, z 2 ) defined in ([7]). (We will also sometimes write G £ (zi,z 2 ) := G £t£ (zi, z 2 ).) 

First we define the function £, z (y), for y E D, to be equal to — logmax(e, \z — y\) minus 
the harmonic extension G z (y) to D of the restriction of — logmax(£:, \z — y\) to dD. Observe 
that this £ £ (y) tends to zero as y — > dD and that as a distribution — A£| (restricted to D) is 
equal to 2irp z , where as before p z £ is a uniform measure on D DdB £ (z). Integrating by parts, 
we immediately have the following: 

Proposition 3.2. The function G eit£2 (z 1 , z 2 ) is equal to the mean value of on the cir- 
cle dB E2 (z 2 ). In particular, if B £l {z\) and B £2 (z 2 ) are disjoint and contained in D then 
G £lj£2 (zi, z 2 ) = G(zi,z 2 ). If B £l (z) C D and Si > e 2 then 

G £1 ,e 2 ( z , z ) = -logei + \ogC(z;D). 

Proof of Proposition I3.ll We first claim that for each e an d D there exists a constant 
K such that 

Var (h ei (zi) - h £2 (z 2 )) < K\e x - e 2 \ + K\z\ - z 2 \ 
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for all zi,Z2 G D and 61,62 G [eo,oo). Since the variance can only increase if D is replaced 
with a larger domain, it suffices to show this holds when D is replaced by a sufficiently large 
disc D' (say, centered in D with 10 times the diameter r of D), and e is restricted to values 
in [eq, 5r]. (For larger values of e, the set dB £ (z) cannot intersect -D when 2; G -D.) Since 

Var {h £l (zi) - h £2 (z 2 )) = G ei>ei (zi, Zi) - 2G £l)£2 (z 1 , z 2 ) + G £2 . £2 {z 2 , z 2 ), 

it suffices to show that G Eue2 (zi, z 2 ) is a Lipschitz function of (ei, e 2 , Z\, z 2 ) for the range of 
(6i,6 2 , zi, z 2 ) values indicated above. This follows from Proposition 13.21 and the fact (whose 
proof we leave to the reader) that £| is a Lipschitz function when z G D and e > Eq, with a 
Lipschitz constant that holds uniformly over these e and z values. 
The claim implies that for some K we have 

E[\h ei (z 1 )-h £2 (z 2 )\ a ]<K\z 1 -z 2 \ a / 2 

for all a > 0. This puts us in the setting of the multiparameter Kolmogorov-Centsov theorem 
[KS91[ IPM] . which states the following: Suppose that the random field X(a), a G niLifO)^] 
satisfies E[|X(a) — X(6)| a ] < K\a — b\ 1+l3 for all a, b, for some fixed constants a, (3, K. Then 
there exists an almost surely continuous modification of the random field and this process is 
7-Holder continuous for every 7 < (3 /a. Applying this for large a allows us to deduce that 
h £ (z), as a function of e and z, is locally 77-Holder continuous for all 77 < 1/2. □ 

Proposition 3.3. Write V t = h e -t(z), and = inf{t : B e -t(z) C D}. If z G D is fixed, 
then the law of 

Vt ■= Hg+t — Hg 

is a standard Brownian motion in t. 

Proof. Since we already know that the h £ (z) are jointly Gaussian random variables, it is 
enough to compute the variances of h £ (z) and h £ >(z) for fixed 6,6', and these are given in 
Proposition 13.21 □ 

3.2 Random metrics: Liouville quantum gravity 

The remainder of the paper makes frequent use of the following simple fact, which the reader 
may recall (or verify): if N is a Gaussian random variable with mean a and variance b then 

Ee N = e a+b/2 . (8) 

Since Kh £ (z) = when h is an instance of the GFF with zero boundary conditions, we have 

Ee-r^W = e Var[ 7 M2)]/2_ 

Recall G £ (z, z) = log C{z; D) — \og6 when B e (z) C D. Then we have 
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More general moments of the random variables e lhe ^ are also easy to calculate. For 
example, we have 

Ee 7My) e 7M*) = enp(VM[y(h E (y)+h £ (z))}/2) = exp f^- (G e (y, y) + G £ (z, z) + 2G £ (y, z)) 

(10) 

By Proposition 13.21 we have G e (y, z) = G(y, z) whenever \y — z\ > 2e and B £ (y)U B e (z) C 
D. In this case we have 

Ee jh e (y) e ihe(z) = ^ C(y;D)C(z;D) y /2 ^ G{y ^y 

Write h £ = yh e + ^- loge. Then we have 

Ee Uz) = C(z;Dy 2 / 2 x 1 

and when \y — z\ > 2e we have 

Ee h s (y) e Uz) = ( C (y; D) C(z; D)f /2 e ^ GM x (C(y; D) C(z; D)f 2 \y - z\~^ x \y - z\~^ 

where x indicates that equality holds up to a constant factor when y and z are restricted to 
any compact subset of D. 

Now, for each fixed e, write fx £ := e he ^dz (which in essence corresponds to the "Wick 
normal ordering" of the original measure [Sim74j ) . We now argue that these converge weakly 
to a limiting random measure on D. 

Proof of Proposition 11.11 Fix 7 G [0, 2). It is easy to see that if for each diadic square 
S compactly supported in D the random variables fi 2 -k(S) converge to a finite limit as 
k — > 00, almost surely, then fi 2 -k almost surely converges weakly to a limiting measure. We 
will prove convergence of fi 2 - k (S) by showing that the expectation of |/j, 2 - fc (S 1 ) — fi 2 -k-i(S)\ 
decays exponentially in k. Without loss of generality, we may assume S is the unit square 
[0, l] 2 , so that He(S) is precisely the mean value of e he ^ on S. 
As shown above, we have 

Ee he(*) = c{ Z] Dyi 2 '\ 

(which is bounded between positive constants) when z G S and e is sufficiently small. 

For y = (yi,y 2 ) £ (0, l) 2 and k > 1, let be the set of 2 2k points (a, b) G S with the 
property that {2 k a — y±, 2 k b — y 2 ) G Z 2 . Let A y k be the mean value of exp h 2 -k-i(z) on the 
set S%. Let B\ be the mean value of exp h 2 -h-%{z) over the same set. 

Clearly, ^ 2 -k-x{S) is the mean value of A v k over y G [0, l] 2 and fi 2 -k-2(S) the mean value 
of B v k over y G [0, l] 2 . Applying Jensen's inequality to the convex function | • |, it now suffices 
for us to prove that E|^4| — B y k \ tends to zero exponentially in k (uniformly in y). Since the 
balls of radius 2~ k ~ 1 centered at points in S| do not overlap, we have that conditioned on 
the values of h 2 -k-i{z) for z G Sf , the random variables h 2 -k-2{z) for z G S k are independent 
of one another; each is a Gaussian of variance log 2 and mean h 2 -k-i(z). 

Hence, given the values of h 2 -k-i(z) for z G S^, the value of conditional variance of 
A y k — B v k is given by 

Var (V^- fe - l(z) - e^- k - 2{z) \h 2 -k-i{z)^ = 2^ k C (eV*- 1 ^) 2 , (11) 
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where 

C = Var (e^'-iW - e h ^ k ^ {z) \h 2 - k -i{z) = o) . 

Note that C is a constant that does not depend on k and z. The unconditional variance of 
A k — B\ is given by the expectation of (TTTT) . It is tempting to argue that this expectation 
tends to zero exponentially in k (which would in turn imply that E\A y k — B\\ tends to zero 
exponentially in k) , but this turns out to be true only for < 7 2 < 2 and not for 2 < 7 2 < 4. 
To see this, note that 



E 



r 2 e h s ( z)) 2 



T 2 loge i -.2 i 



Summing over the e -2 points yields (up to an e-independent constant factor) e 2 ^ 2 , which 
does not tend to zero when r y 2 > 2. 

However, we can deal with the case 7 2 > 2 by breaking the sum defining A y k — B v k into two 
parts and dealing with them separately. The idea is that the variance of A y k — B y k is only large 
because of rare occurrences (where h e (z) is much larger than typical) whose contribution to 
the expectation of A v k is exponentially small in k. 

To make this precise, fix some a with 7 < a < 27 and let Sf denote set of points z E S k 
with the property that h £ (z) > — aloge, where e = 2~ k . Let A y k denote the mean value of 
Igv exp h 2 -k-i(z) and B y k the mean value of lgv exp h 2 -k-2(z). 

We claim that EA y tends to zero exponentially in k. To see this, observe that for fixed 
z £ S, the random variable h e (z) is a centered Gaussian with variance a 2 = — loge plus a 

constant; and the expectation of e he ^ — which we know to be constant for all e small enough 

_ 1 2 

so that z is distance at least e from the boundary of D — takes the form J R e z^e^drj. This 
expression is a constant times a Gaussian integral of mean 7a 2 and variance a 2 ; thus 

is the probability that a normal random variable 77 of mean 70" 2 and variance a 2 satisfies 
i] > aa 2 , and this clearly tends to zero exponentially in a 2 , from which the claim easily 
follows. 

Note that by Proposition 13.31 E5^ = EA^, and in particular Ei^ also tends to zero 
exponentially in k. Since B y k and A y k are non-negative, applying Jensen's inequality to | • | 
shows that E|i^ — A y k \ tends to zero exponentially in k. 

For the next step, we wish to bound E|(_B^ — B k ) — (A k — A y )\, which requires us to 
consider points in S k \ S k and the expectation 



e 4 E 



I , I P h e {z) 

± r/<acr z I 



e'Eil^e 2 ^ 2 ^]. (13) 



This differs from (j!2p by a factor that represents the probability that a Gaussian random 
variable with variance — loge and mean —27 loge is less than —aloge. Since a < 27, this 
probability decays exponentially in — loge at rate (27 — a) 2 /2. Thus (jTSJ) becomes, up to a 
constant factor (universal in e and z £ S), 



2 £ (2 7 -a) 2 /2_ 
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Summing over the e 2 points, we obtain 

£ 2- 7 2 + (2 7 -a) 2 /2_ 

To conclude, we only need to make sure we chose a G (7, 27) small enough so that the sum 
in the exponent is positive, and this is clearly possible. In fact, taking a close to 7, the 
exponent becomes close to 2 — 7 2 + 7 2 /2 = 2 — which is positive when 7 < 2. □ 

3.3 Rooted random metrics 

Before proving Proposition ll.2[ we introduce a notion of rooted random metric and use it to 
prove a uniform integrability result for the random variables fi £ (S) discussed above. 

Write £ := Z £ e^ {z) dzdh, where Z £ is a constant chosen to make £ a probability 
measure. Sampling from e may be described as a two step procedure. First sample z from 
its marginal distribution. Then sample h from the distribution of the Gaussian free field 
weighted by e 7/l£ ^- ) . The latter has the law of the original GFF plus 7^ where satisfies a 
Dirichlet problem: — A£* is the multiple of the uniform measure on dB £ (z) with total mass 
2tt (because h is yht times the standard GFF; if h were the standard GFF the total mass 
would be 1). As noted in Section I3TT1 this £| has been computed explicitly: 

C £ (y) = -logmax{\z - y\,e} - G K (y), 

where G z is the harmonic interpolation to D of the first term on dD, as long as B e (z) C D. 

For each fixed e, the marginal distribution of z is given by f(z)dz where f(z) = Khe lh ^ z '. 
Thus f(z) is proportional to e^ / 2 ) 1o sC(z-,d) by q 

Let be the limit of the measures e as e — > 0: that is, is the measure on pairs (z, h) 
for which the marginal distribution of z is proportional to e^ 7 i 2 > lo ^ c \ z > D >dz and, given z, the 
conditional law of h is that of the original GFF plus the deterministic function 7^ (viewed 
as a distribution). When S C D we will also write Q £ for the measure £ conditioned on 
the event z G S. The following is obvious from our definitions: 

Proposition 3.4. With probability one, z is a 7-thick point of h. That is, 

liminf h £ (z) /\oge~ l > 7. 

e— >0 

(In fact, the limit exists and equality holds almost surely.) 

Since the marginal law of h is absolutely continuous with respect to the law of h (with 
Radon-Nikodym derivative this implies that \ih is almost surely supported on 7-thick 

points. It was shown by Hu, Miller, and Peres that the set of 7-thick points has Hausdorff 
dimension 2 — ^ almost surely jHMPj . 

Proof of Proposition 11.21 The almost sure weak convergence of the fi n to a limit jl is 
immediate from the martingale convergence theorem. Recall the expression ([T]) 

H n = exp (-yh n (z) - ^-Varh n (z) + \ogC(z; D)\ dz, 
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and observe that for each z, the exponential term 



(2 2 
jh n {z) - —Vaih n {z) + — logC(z; D] 
2 2 

is a non- negative martingale with respect to the filtration of h n . (This is a consequence 
of (JE])-) Fubini's theorem implies that [i n {A) is a martingale for any Borel measurable set 
A C D, and the martingale convergence theorem implies that the limit lim fi n (A) exists 
almost surely. In particular, this holds whenever A is a diadic square contained in D and 
from this easily follows the desired weak convergence. 

We still need to show that jl = fi almost surely, where /i is as constructed in Proposition 
11.11 Let h™ denote the mean value of h n on dB £ (z). For each particular choice of z, and e 
small enough so that B £ (z) C D, and for each n, we have 

(2 2 \ 

7/£(z) - yVar^(z) + -L log C(z; D) J . 

Taking the limit as e — > and using the continuity of and Var h n (z) and the expression 
([I]), we have 

limE[// e (S)|fc n ] =// n (5) (14) 

for each diadic square S. 

We next wish to argue that 

limE[fi £ {S)\h n ] = E[lim fi £ {S)\h n ] = fi n {S), (15) 

provided that < 7 < 2. 

We first argue this in the case n — 0. Since Proposition 11.11 implies the existence of the 
limit, it is enough to show that the random variables M e = ^ E (S) are uniformly integrable as 
e — > 0. Let M = EM e for e small enough so that B e (S) C D. (By (jHJ) this expectation is the 
same for all sufficiently small e.) The uniform integrability is equivalent to the statement 
that the probability measures r] e := M~ l M e dM e are tight, i.e., for all 5 there exists a constant 
C > such that r) £ ([C, 00)) < S for all e. (Here M~ 1 M £ dM E denotes the probability measure 
on R whose Radon-Nikodym derivative with respect to the law of M £ is given by M~ 1 M £ .) 
Since M £ is a function of h, this is equivalent to the statement that with respect to the 
measure M~ l M £ (h)dh the random variables M £ (h) are tight. This in turn can be rewritten 
as the statement that for each S we can find a C such that 

Q s £ {M £ {h) >C}<5 

for all e. 

Let e = sup-jV : B £ >(S) C D}. We sample the pair (z, h) in the following steps. In the 
first step, we sample z from its marginal law (which is independent of e for e sufficiently 
small). Write h = h — 7^*. Given z, the law of h is that of a GFF on D. In the second step, 
we sample B t = h e -t £o (z) — h £Q (z) for all t G [0, — log(e/eo)]- By Proposition 13. 3^ B t is a 
Brownian motion on this interval independent of z. The conditional expectation of h given 
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the whole process B t (which we have defined only for t £ [0, — log(e/e )]) and z is given by 
the function (viewed as a distribution) 



[ B , max.{\z-w\,e\ \Z — W\ < Bn 

ti(w):=E[h(w)\z,Bt] = \ ~ log ^ . 

[ | z — W | > Sq 

Note that once z is fixed, for each w the mean value of h'(-) on dB £ (w) (which we 
denote by h! £ {w)) is a weighted average of B t over values of t between V — log ^~^ +£ and 

— log ^~^~ £ A — log^. From this it is not hard to see that given z the variance of h' £ (w) 
is between these two values of t. We claim that these bounds differ from each other by at 
most the additive constant log 3. Exponentiating the bounds and multiplying by Eq, this 
is equivalent to the statement that e A (\w — z\ + e) differs from (\w — z\ — e) \/ e by & 
multiplicative factor of at most three (recall that e < e ), which is easily checked. Thus the 
variance of h' £ (w) is within log 3 of the former bound 



u = u(w) = V — log 



\W — Z\ + E 



The conditional variance of h £ (w), given h' £ (w), is within log 3 of its initial variance minus 
this value. Thus, with respect to 0f , we have 

E[ £ T 2 /2 e 7M«0|^ St] x exp (~i~ti £ {w) + 7 2 ^(w) - 7 2 m(w)/2) x exp (^~h' £ (w) + -f 2 u{w) /2^j , 

where x indicates equality up to a multiplicative factor bounded between positive constants 
uniformly in e and z. 

Now, given any positive constants a and b, there is a positive probability that a Brownian 
motion B t run for an infinite amount of time will satisfy ^B t < a + bt for all t > 0. In fact, for 
each fixed 6, this probability can be made as close to one as possible by taking a sufficiently 
large. Since < 7 < 2 we can choose a value of b with < b < 2 — 7 2 /2. Then note that 
when ^Bt < a + bt for all t we have, for some constant C 

E[E^ 2/2 e^ w) \z,B t } < C e a \z-w\' b ^ 2/2 , 

for \z — w\ < £0, which in turn implies that 

E[n £ (S)\z,B t ] < [ C G e a \z-w\- h -^l 2 dw + E\}i £ {S\B £0 {z))] 

JB eo {z) 

< [ C Q e a \z-w\^ 2/2 dw + E[^ £ {S)l 

JB eo (z) 

and since b + 7 2 /2 < 2 (and E[fi £ (S)] is constant for sufficiently small e), the right hand 
side is at most a finite constant C\ = C\{a) that is independent of e. Now, given b and a 
constant 5 > we can choose a large enough so that the probability that 7^ < a + bt for 
all t is at least 1 — 5/2. Then we take C = c ^ . If there were probability at least 5 that 
/i £ (S) > C then there would have to be probability at least 5/2 that jB t < a + bt for all t 
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and fi £ (S) > C, which could contradict our bound on the conditional expectation of fi £ (S) 
given that jBt < a + bt for all t. This implies that the probability measures t] £ are tight, 
which in turn completes the proof of ( fT5l) in the case n = and h° = 0. 

To extend to the case n ^ 0, note that since the random variables /^ £ (S) almost surely 
converge to a limit (with expectation lim^o Efi £ (S)), it must be the case that conditioned 
on h n (for almost all h n ), we still have that fi £ (S) almost surely converges to a limit; the fact 
that 

Epim Li £ (S)\h n ] < limE[/z £ (S)|/i n ] 

for almost all h n is immediate from Fatou's lemma, and the fact that equality holds almost 
surely follows from the unconditional result. The extension of (fl5|) to non-zero h° is trivial 
for functions that are piecewise constant on diadic squares, and the more general case follows 
easily by approximation by piecewise constant functions. 

Proposition 11.21 is an immediate consequence of (THI) and (TTBT) . □ 



4 KPZ proofs 

4.1 Circle average KPZ 

Fix z G D and some Sq such that B eo (z) C D. For any e < eq write t = — \og(e/eo) and 
Vt = h E (z) — h £0 (z). The law of Vt is that of a Brownian motion with Vq = (by Proposition 
13. 3p . It follows from Proposition 11.21 and recalling the notation of Proposition ll.il that the 
expectation 



E h [fi h (B £ (z))\h £ (z)]=E h 
has approximately the form 



e^dzlVt 



in the sense that the ratio of the logarithms of the two quantities tends to 1 as e — > 0. 

Definition 4.1. Let B s (z) be the largest Euclidean ball in D centered at z for which (To)) is 
equal to 5. The radius of this ball is e~ TA where 

T A := inf{t : -V t + Qt = A}, 

and A := —(\og5)/j. 

As a step towards Theorem 11.51 we prove the following in this section, which is perhaps 
the most straightforward form of KPZ to prove: 

Theorem 4.2. Theorem \ 1.51 holds with B s (z) replaced with B s (z). That is, in the setting of 
Theorem \1.5l if 

\ogEfi {z : B £ (z) e X} _ 
lim — x, 

e^o log e 2 
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then it follows that, when X and fi are chosen independently, we have 

\ogEn{z :B 5 (z) G X} 

lrm — = A, 

5^0 log 5 

where A is the non-negative solution to 



We present two proofs: the first based on exponential martingales, the second based 
on large deviations theory and Schilder's theorem. (The first proof is shorter, but readers 
familiar with large deviations of Brownian motion will recognize that it is essentially the 
second proof in the disguise.) 

Both proofs use the fact that 

E h fx{z : B\z) G X} 

is proportional to 

Q{(z,h) : B\z) G X}, 

to replace an expectation computation with a probability computation. (Recall the definition 
of O from Section 1X51 ) While this rephrasing is not strictly necessary for the expectation 
computation below, it is conceptually quite natural. 

We use the definition of Vt given above, and assume that the fixed Eq is smaller than the 
distance from D (recall that this was the compact subset of D in Theorem I1.5P to dD. 

As mentioned in Section 13.31 the G conditional law of h given z G D is that of the 
original GFF plus the deterministic function —7 log \z — y\. Thus (for z restricted to points 
of distance at least £0 from dD) the conditional law of V t given z is that of B t + jt, where 
Bt evolves as a standard Brownian motion — in particular, z is independent of the process Vt. 

Proof. The O law of Ta is that of 

inf{t : B t + at = A = -(log<5)/ 7 }, a := Q - 7 = - - ? > 0, (17) 

7 2 

where (±)£>t is standard Brownian motion with Bq = 0. Since z is independent of Ta, the 
theorem hypothesis implies that conditioned on Ta, the probability qA that the ball of radius 
e~ TA centered at z is in X is approximately exp (— 2xTa), in the sense that the ratio of the 
logs of these two quantities tends to 1 as Ta — > 00. Computing the expectation 

E[exp(-2xT A )], (18) 

with respect to a random T4 will give us upper and lower bounds on qA since it easily follows 
that 

E [exp (-2x 1 T A )} <q A <^ [exp (-2x 2 T A )} , (19) 

for any fixed < x<i < x < x\ and sufficiently large A. 

To compute (TT81) . consider for any (3 the exponential martingale exp(j3B t — f3 2 t/2). At 
the stopping time Ta 

E[eMPBT A -P 2 T A /2)] =1. 
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By definition Bt a = A — clTa- Thus, 

Eexp[-(pa + P 2 /2)T A ] = exp(-/M). 
Setting 2x := Pa + /3 2 /2, we obtain 

Eexp(-2a;T A ) = exp(-/M). 



(20) 



Now if we set A = fi/j, and a = Q — j = - — |, we find that the equation 2x := (3a + (3 2 /2 
is equivalent to the KPZ formula. The continuity of this expression and ffl9|) together yield 
the theorem. □ 
We remark that the above yields the explicit probability distribution Pa(£). The inverse 
Laplace transform Pa(£) of := Eexp(— 2xTa), with respect to 2x, is the probability 

density such that PA{t)dt := Prob (Ta £ [t,t + Its explicit expression is 



P A (t) = (27r)- 1/2 Ar 3/2 exp 
where as above we have A = — (log5)/7, t = ■ 



-(1/2) {Ar 1 ' 2 -at 1 ' 2 )' 
log £ and a = Q — 7. 



(21) 



4.2 Large deviations proof of circle average KPZ 



In this section, we present an alternative proof of Theorem I4.2[ using Schilder's theorem. 

Lemma 4.3. Fix a constant a > 0. Let B t be a standard Brownian motion. For each A > 0, 
write 

T A = inf{t : B t + at = A}. (22) 

Then the family of random variables A~ x Ta satisfies a large deviations principle with speed 
A and rate function 

I(rj) = v (^--a^ j2 = r j - 1 /2 - a + a 2 rj/2. 

Proof. Schilder's Theorem (see Theorem 5.3.2 of [DZ] ) gives an LDP for the sample path 
of a~ x Bt (where Bt is standard Brownian motion) with speed a 2 and rate function given 
by the Dirichlet energy. The variable A~ x Ta can be written as inf{t : Wt + at = 1} where 
W t = BAt/A, which has the same law as V ' A~ x B t . Clearly, among all functions <\> G #i([0, 00)) 
satisfying 0(0) = and inf{t : <p(t) + at = 1} < 77, the one with minimal Dirichlet energy is 
given by 

[\-a)t t<r] 
[~-a)v t>V- 

By the contraction principle (Theorem 4.2.1 of [DZ] ). the rate function desired in Lemma 
14.31 is given by this minimal Dirichlet energy, i.e., I{rf) = rj(- — a) 2 /2. □ 

Lemma 4.4. Consider the following two part experiment. First choose Ta as above. Then 
toss a coin that comes up heads with probability 

e~ 2xT \ 



4>{t) 
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Then the probability that the coin comes up heads decays exponentially in A at rate (3 where 
(3 and x are related by 

(3 = inf {I(rj) + 2xr]} , (23) 
v 

or equivalently by 

Ax = f3 2 + 2a/3. (24) 

Proof. The exponential decay with the exponent given in (1231) is an immediate consequence 
of Varadhan's integral lemma (Theorem 4.3.1 of |DZJ). To derive (I24p from (123ft . we set the 
derivative of I(rj) + 2xrj to zero and find —r]~ 2 /2 + a 2 /2 + 2x = 0. Hence the minimum is 
achieved at 

r] = (a 2 + Ax)- 1/2 . (25) 
We then compute (3 = I(r)o) + 2xrj to be 

(a 2 + Ax) 1/2 /2 - a + a 2 (a 2 + Ax)- 1/2 /2 + 2x(a 2 + Ax)- 1 ' 2 . 

Simplifying, we have (3 = (a 2 + Ax) 1 / 2 — a, which is equivalent to (I24j) . □ 

Proof of Theorem 14.21 As above, we aim to show that P{B s (z) E X} scales as 
e _/3A = = 5 A where A = P/j, where 5 and e are related via the stopping time Ta (|17p . 
Rescaling Ta by A" 1 as in (1221) puts us in the framework of large deviations Lemma \A. 31 As 
above, to describe the probability P{B 5 (z) G X} we can imagine that we first choose the 
radius e of B s (z) and then toss a coin that comes up heads with probability e 2x to decide 
whether the ball is in X . This puts us in the framework of the second large deviations 
Lemma [4.41 Using (I24p . we have 

Ax = (3 2 + 2a(3 = (7A) 2 + 2a 7 A, 

where a = Q — 7. Plugging in this value of a and simplifying, we obtain the KPZ relation 

x = \ ( 7 2 A 2 + 2 7 (Q - 7 ) A) = ^A 2 + (l - A- 

As in the previous proof, if the probability given e is not exactly e 2x , but the ratio of the 
log of this probability to the log of e 2x tends to 1 as e — * 0, we obtain the same theorem by 
using alternate values of x to give upper and lower bounds. □ 

The optimum r] = (a 2 + Ax)~ x l 2 obtained in (125!) has a natural interpretation — it 
suggests that (in the large deviations sense described above) T A /A is concentrated near r] . 

Equivalently, since A = - = ^ +4a ^ 1 ~ a , we can say that A/Ta is concentrated near 

7 A + a = 7 A + Q — 7, which implies that is concentrated near 7 ( 7 A + Q — 7). Note 

lOg £ 

that the same result can also be obtained directly from the explicit probability density (121]) . 
This is the concentration one obtains at an a-thick point of the GFF h, where 

a = 7 - 7A. (26) 

Very informally, this suggests the quantum support of a quantum fractal of dimension A is 
made up of a-thick points of h. This generalizes the idea of Proposition 13 A\ which concerns 
the case A = 0. 
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4.3 Probability total mass is very small 

Lemma 4.5. LetH) = Pi(0) be the unit disc and fix 7 G [0, 2) and take [i = e lh ^dz as defined 
previously. Then the random variable A = log fi (-81/2(0)) satisfies Pa{v) := < i]] < e~ c ' v 
for some fixed constant C > and all sufficiently negative values off]. 

Proof. Let h' be the projection of h onto the space of functions in if (O) that are harmonic 
inside the two discs Pi/ 4 (l/4) and i? 1 / 4 (— 1/4). Recall that the orthogonal complement of 
this space is the space of functions supported on these discs, or more precisely, the space 
if [.61/4(1/4) UPi/ 4 (— 1/4)]. Hence, the law of h — h' is that of a sum of independent Gaussian 
free fields on Pl/ 4 (1/4) and Pi/ 4 (— 1/4) with zero boundary conditions. 

Let h be the infimum of h! over the union of the two smaller discs P_ = Pi/g(— 1/4) and 
B+ = Bi/ 8 (l/4). Write A_ = log/i ft _ ft >(P_) and A + = Iog^ h _ h /(S + ). By Proposition O the 
law of each of A + and A_ is the same as the law of A + 7Q log(l/4) = A — 7Q log 4; clearly 
A + and A_ are independent of one another. Also, fi{B + ) > e~ t -fif l -h'{B+) (and similarly for 
B_), which implies 

A > max{A„, A + } + jh. (27) 

First we will show that the probability distribution of h has superexponential decay. Since 
h! is harmonic on B + (with ft/ (1/4) = hy 8 (l/4)) this h' is the real part of an analytic function 
on B + . In particular, h' restricted to B + can be expanded as /i'(l/4) + Y^=i R- e [ a n^ n (z — 
1/4) n ] for some complex a n . Since each of the random variables Rea n and lma n is a real- 
valued linear functional of h', it is a Gaussian random variable. 

Let F be any linear functional on subspace H'(D) C H(D) of functions that are harmonic 
on B + ; then F(h) can be written as (h, /)y for some / G H'(D) and the variance of -F(ft) 
is (/, /)v- Note that (/, /) v x is the smallest Dirichlet energy obtained by F on the set of 
functions {g G H'(D) : -F(g) = 1} (since / has minimal energy on the set {g G H\D) : 
/••(.'/) (/-./>}). 

In the case of the linear functional Rea n (the case lma n is similar), this minimal Dirichlet 
energy is at least the Dirichlet energy of Re4 n (z — l/4) n restricted to i?i/ 4 (l/4). (If for any 
g G H(D), the linear functionals Ima n or Rea m for m 7^ n applied to g are non-zero, then 
the Dirichlet energy of g restricted to i?i/ 4 (l/4) will be greater than if they were zero, by 
orthogonality of Re z n and Im z m on D; and the total Dirichlet energy of g on D is at least 
the Dirichlet energy of the restriction to £?i/ 4 (l/4).) By conformal invariance of the Dirichlet 
inner product, this energy is given by 

/ ' n 2 \z n ~ 1 \ 2 dz = n 2 [ r 2n - 2 2nrdr = nn. 
Jb Jo 

This implies that the variances of Re a n and Im a n are at most 1 / (nn) . 

In particular, the variance of |a„|r n , for any fixed r < 1, will decay exponentially in n. 
Thus, the probability that even one of the a n satisfies |a„|r ra > c, where c is a fixed constant, 
decays quadratic-exponentially in c. It follows that the probability density function p of h 
satisfies p(rf) < e~— v2 for some C_ > and all sufficiently negative 77. 

Now, let Pi (77) be the probability that h < .I77/7 and A < 77. Let P2(v) be the probability 
that A < 1] and h > .177/7. Then Pa{v) = ^[A < 77] = Pi + P2. From the above discussion, 
we have Pi (77) < e~— n for all sufficiently negative values of 77. Note from fl27j) that 
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P2(v) < [pA(.9r/ + 7 gio g 4)] 2 

and 

PM < [Pi(.9r ] + 1 Q\ogA) + P 2 (.9r ] + 1 Q\ogA)] 2 < [e" cV + P 2 (.9r] + 7 Q log 4 

for some C Fix a sufficiently negative 770 and inductively determine r/k via r/k-i = .9% + 
7<5 log 4. The above can be stated (with a modified C) as 



p 2 m< (e- c '^-+p 2 (%-i; 



is 



If we write p*, = <i , then this can be restated as pk < (1 + £>fc-i) 2 e 2C ( -' 7fc - 1 ^Z 2 ). It i 
easy to see that we can have pk > 1 for only finitely many k, which implies that the lemma 
holds when restricted to the sequence Because of the monotonicity of pa(v)i this implies 
the lemma for all 77. □ 
From the lemma above, it is easy to derive the following, which includes a restatement 
of (fTEl) . together with a strong upper bound on the probability that fi(B £ (z)) is much lower 
than this expectation. 

Lemma 4.6. Fix z and e so that B £ (z) C D. Then 

E^(B £ (z))\h £ (z)]^s^e^\ 

where 

2 7 
7 2 

as in Proposition \2.1\ Moreover, conditioned on h £ /(z), for all e' > e, we have that 



P 



^B £ {z)) e _ t 

£ lQ e lh s {z) 



< Cie 



-c 2 t 2 



for some positive constants C\ and C 2 independent oft>0, z, D, and the values h £ >(z) for 
e'>e. 

Roughly speaking, the above lemma says that the total quantum mass in a ball is unlikely 
to be a lot smaller than the mass we would predict given the average value of h on the 
boundary of that ball; the following says that (even when we use the measure), the total 
quantum mass has some constant probability to be (at least a little bit) smaller than this 
prediction. 

Lemma 4.7. Let z and h be chosen from s for a fixed compact subset S of D, and fix a 
5 > 0, with quantum balls B 5 (z) and B 5 (z) defined as in Definition \1.3\ and Definition 



Conditioned on the event B (z) C D and on the radius of B (z), the conditional probability 
that B 5 {z) C B 5 (z) is bounded below by a positive constant independent of D, D, and S. 
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4.4 Proof of interior KPZ 



In this section we derive Theorem 11.51 as a consequence of Lemma 14.61 Lemma 14.71 and the 
arguments in Theorem 14.21 

Proof of Theorem II. 5L We use the same notation as in Theorem 14.21 but we write 
T a = — loge where e is the radius of B s (z). The proof of the Theorem 14.21 carries through 
exactly once we show that 

logEjexpM^)] = 
a^oo logE[exp(-2xT j4 )] v ; 



Lemma I4.6| applied for t/7 = A(l — a) and for a fixed a < 1, implies that given Ta, the 
probability that T a < T a A decays superexponentially in A for any a < 1. This implies that 

logE[exp(-2xT A )] < i 

a-,00 logE [exp (-2xT aj4 )] ~~ 
and since this holds for all a < 1, the result follows immediately from Lemma [4.71 and the 
continuity of the coefficient of A in the exponent in ( 1201) . □ 



5 Box formulation of KPZ 

In this section we prove Proposition 11.61 

Proof of Proposition 11.61 The first fact is standard; observe that if e is a power of 
two then S £ (X) C B 2e (X), since the ball of radius 2e about a point contains any diadic 
box of width e that contains the same point. Similarly, B 2e (z) is contained in the union of 
a diadic box — of width 2e, containing z — with the eight diadic boxes of the same size 
whose boundaries touch its boundary. This implies that B 2e {X) is contained in the union of 
S 2e (X) and corresponding 8 translations of S 2£ (X), so fj,o(B 2E (X)) < 9fio(S 2e (X)). 

For the second part, we first prove the statement where we replace 5N(n, 5, X) with the 
total mass of the set S S (X) of points in the diadic boxes of quantum area at most 5 that 
intersect X. In this case, the proof is similar to that of Theorem 14.21 We use the notation 
of Theorem 14.21 but set Ta to be — log£ where e is the largest value for which the diadic box 
with edge length e has \i area at most 5. The remainder of the argument is essentially the 
same as the proof of Theorem 11.51 We need only replace Ta with Ta and note that 



log 


E 


exp 


-2xf A ) 


lo£ 


;E 


exp 


-2xT A )\ 



= 1, 



by the same argument used to compare Ta and Ta- 

Having obtained this result, it remains only to show that 

logE<$7VQM,X) 
™ logE/^pO) 

The liminf is clearly at most 1, since by definition, the denominator is greater or equal to the 
numerator. The other direction uses the same argument at the end of the proof of Theorem 
11.51 It suffices to note that if we fix a constant a < 1, Lemma 14.61 implies that given Ta, the 
probability that Ta < T a A decays superexponentially in A for any fixed a < 1. □ 
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6 Boundary KPZ 



Most of the results in this paper about random measures on D have straightforward analogs 
about random measures on 3D. The proofs are essentially identical, but we will sketch the 
differences in the arguments here. 

Suppose that D is a domain with piecewise linear boundary and that h is an instance of 
the GFF on D with free boundary conditions, normalized to have mean zero. 

This means that h = Yl n a nfn where the a n are i.i.d. zero mean unit variance normal 
random variables and the /„ are an orthonormal basis, with respect to the inner product 

(/ 1 ,/ 2 ) v :=(27r)- 1 / Vf 1 (z)-Vf 2 (z)dz, 

J D 

of the Hilbert space closure H(D) of the space of C°° bounded real- valued (but not necessarily 
compactly supported) functions on D with mean zero. 

Note that if / is a compactly supported smooth function on D for which —Af = p, then 
integration by parts implies that the variance of (h, p) is the Dirichlet energy of / — same 
as in the zero boundary case. Similarly, suppose that / is a smooth function that is not 
compactly supported but has a gradient that vanishes in the normal direction to 3D, and 
we write p = —Af. Then integration by parts implies that the variance of (h,p) is (/, /)y- 

We can also make sense of h £ (z), for a point z on the boundary of D, to be the mean value 
of z on the semicircle of radius e centered at z and contained in the domain D. In this case, if 
Eq is small enough so that B £o contains no corners of D and exactly one semi-disc of B £o lies 
in D, we have that h £ (z) — h £o (z) is a standard Brownian motion in time 2t = — 21og(e/eo)j 
as in Proposition 13.31 The 2t in place of t comes from two factors: first, by integration by 
parts, h e (z) — h eo (z) is equal to (h, where £(•) is the continuous function which is equal 
to 2 log \z — -| on the half-annulus HI D {y : e < \y — z\ < 6q} and is constant outside of the 
half-annulus. (The 2 log \z — -| in place of log \z — -| comes from the fact that we are taking 
an average over half a circle.) The variance of h £ (z) — h £0 (z) is given by the Dirichlet energy 
(£>0v5 which is twice as big as before (given the factor of 2 in the definition of £ and the 
fact that the integral is only over half as much area). 

Thus at a point on the interior of one of the boundary lines of D, the variance of h £ (z) 
scales like —2 loge instead of — logs. We define the boundary measure pf := e 1 ^ 4: e yh ^ z ^ 2 dz, 
where in this case dz is Lebesgue measure on the boundary of D. Here we use e lhe ^ z ^ 2 instead 
of e 7he ^ because we are integrating a length instead of an area; as before, the power of e 
chosen makes the factor preceding dz an exponential martingale in time t = — loge. 

We define p B to be the weak limit as e — > of the measures pf (see the theorem below 
for existence of this limit when < 7 < 2). For z G 3D we write B £ (z) := B £ (z) H dD and 
we define B s (z) to be the (largest) set B e (z) whose p B measure is 5. 

Likewise define 

B £ {X) = {zedD: B e {z) n X ^ 0} 

and 

B\X) = {zedD: B 5 (z) n X ^ 0}. 

We say that a (deterministic or random) fractal subset X of the boundary of D has Eu- 
clidean expectation dimension 1—x and Euclidean scaling exponent x in the bound- 
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ary sense if the expected measure of B e (X) decays like e x , i.e., 



e^O log e 

We say that X has boundary quantum scaling exponent A if when X and fi B 
defined above) are chosen independently we have 



as 



lim 

<5^0 



\ogEfi B (B 6 (X)) 
log 5 



A. 



Theorem 6.1. Given the assumptions above, Proposition ] 1.1\ and Theorems \1.5\ and \4-^\ hold, 

precisely as stated, when /i £ is replaced by fif , /i is replaced by \i B ; fio (Lebesgue measure 
on D) is replaced by Lebesgue measure on 3D; B £ and B s are replaced with B e and B s , 
respectively; and the compact subset of D is replaced with a closed subinterval of one of the 
boundary line segments of D. 

Proof. The proofs in the boundary case proceed exactly the same as in the interior point 
case, up to factors of 2 in various places. We sketch the proof of an analog of Theorem 14.21 
in order to indicate where those factors of 2 appear. 

Write t = — \og(e/e ), and let V t = h £ (z) — h eo (z). It is not hard to see that the 
expectation of the boundary line integral 





~^ B (B £ (z))\h £ (z) 




[ e^ h/2 dz\V t 











has approximately the form (which replaces (JT6l) 

exp 



-V t - -Qt 
2 2 . 



(29) 



in the sense that the ratio of the logs of the two quantities tends to 1. Let B s (z) now be the 
largest Euclidean ball B e (z) in D centered at z G dD for which f[2"9~j) is equal to the quantum 
length S. 

As before, the G conditional law of h given z is that of the original GFF plus the 
deterministic function —7 log \ z — y\ (minus a bounded function of y). 

Then given z G dD, the 6 conditional law of V t is that of B 2 t + lt, where B 2 t evolves as a 
Brownian motion with twice the variance of standard Brownian motion, because of the free 
boundary conditions on 3D. Using (1291) . we have 



^ h,2 dz\V t 



UB £ (z) 



exp 



1 n 

-B 2t 



(30) 



This will be equal to the quantum boundary length 8 at the smallest t for which ^ 2 t + r yB 2 t — 
7Qt = 2 log 5. That is, -B 2t + (Q - 7)* = -2 (log 5) ft. If we set A := -(log 5)/ 7, this 
smallest time is a stopping time Ta such that 

T A = infO : B 2t + at = 2A = -2(log<5)/ 7 }, a = Q - 7 = - - J > 0. (31) 

7 2 
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As above, we consider the two part experiment in which we first sample Ta and then sample 
z and check to see whether the ball of radius e = c~ Ta intersects X on the boundary. Given 
T A , the ratio of the logarithms of this probability and 

E [exp {-xT A )\ 

tends to 1 as A — > oo. Consider next for any (3 the exponential martingale exp (§#2* — 
such that 



E 



exp 



2 4 



At the stopping time Ta in particular: 



E 



exp 



E 



P 2 



exp 







Br 



1. 



B 2Ta - —T 



1. 



2 4 

By definition £>2T A = 2 A — clTa- One thus gets the identity 

Eexp[-(/3a/2 + /3 2 / 4 ) t a] = exp(-/M), 

and it now suffices to identify 2x := /5a + /3 2 /2 to obtain the boundary KPZ with A := f3/j 
and 

Eexp(-xT A ) = 5 A = exp(-/M) = exp {-A[{a 2 + 4£) 1/2 -a]}. 

□ 

The reader may observe that the boundary measures described above are preserved under 
the transformations described in Proposition 12.11 One can use this to define the boundary 
measure on more general domains, which may not have piecewise linear boundary conditions. 

We also remark that a similar procedure to that above allows us to make sense of measure 
restricted to lines in the interior of the domain. 



7 Discrete random surface dimensions and heuristics 

Historically, one of the uses of the KPZ formula has been to make heuristic predictions about 
the scaling exponents of random fractal subsets of the plane (see, e.g., [Dup99b , IDFGG OO, 



Dup04 , DupOO , Dup06| , and the references surveyed therein for much more detail). 



In this subsection, we give a very rough and very brief sketch of what such a heuristic 
might entail in a simple example. Readers familiar with discrete quantum gravity models 
(a.k.a. random planar map models, random quadrangulation models, etc.) should note that 
these models have natural interpretations as continuum random metric spaces as well. For 
example, a random planar quadrangulation M n on the sphere — chosen uniformly from 
the set of all simply connected planar quadrangulations with n quadrilaterals — can be 
viewed as a manifold by endowing each quadrilateral with the metric of a unit square. (Of 
course, the resulting manifold will have singularities: negative curvature point masses at 
vertices where more than four unit squares coincide and positive curvature point masses 
at vertices where fewer than four unit squares coincide.) We may then choose a uniform 
square from among this set. Taking an "infinite volume limit" (as n — > oo) one obtains 
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an infinite random quadrangulation M ra with a distinguished square. (See, e.g., [AS03| for 
a precise description of this construction for triangulations.) This infinite random metric 
can be conformally mapped to the plane in such a way that the center of the distinguished 
square is mapped to the origin and the volume of the image of the distinguished square is 
a constant 5 (with a rotation chosen uniformly at random). The images of the unit squares 
of Mqo form a tiling of C by "conformally distorted" unit squares. Different squares have 
different sizes with respect to the Euclidean metric on the plane; intuitively, one would expect 
such a tiling to look something vaguely like the tilings in Figures [U [2j and [3] except that the 
"squares" would be randomly oriented and distorted. The pullback of the intrinsic metric 
of Mqo to the plane via this map takes the form e x (dx 2 + dy 2 ) for some function random A 
(which has logarithmic singularities at the images of the vertices of the squares). Although 
the equivalence of Liouville quantum gravity and discrete quantum gravity is taken as an 
Ansatz throughout much of the literature, to our knowledge the following is the first precise 
conjecture for the complete scaling limit of a discrete quantum gravity model: 

Conjecture 7.1. As 5 — > 0, the function A converges in law (e.g., w.r.t. to the weak topology 
on the space of distributions on the plane modulo additive constants ) to j(h — 7 log | • |) where 
h is an instance of the whole plane Gaussian free field (defined up to additive constants) and 



We further conjecture that other values of 7 are obtained by choosing a random quad- 
rangulation together with a statistical physical model on the quadrangulation (FK cluster 
model, percolation, O(N) model, uniform spanning tree); in this case, the probability of 
a given quadrangulation is proportional to the partition function of the statistical physics 
model on that quadrangulation. (See the references on random matrix theory and geo- 
metrical models cited in the introduction for much more detail; see |Dup06| for a review 
with additional references.) One can also consider scaling limits on spheres or higher genus 
surfaces, as well as different kinds of marked points (corresponding to different logarithmic 
singularities in the scaling limit); however, these are a bit more complicated to describe, so 
we limit attention to the infinite volume case for now. 

By the usual conformal invariance Ansatz, it is natural to expect that if one conditions 
on the infinite quadrangulation, and then samples the loops or trees in these models (as 
mapped into the plane), their law (in the scaling limit) will be independent of the metric. 

Now suppose that for each n we define a random subset X n of M n (for example, X n could 
be the set of the squares hit by a simple random walk started at the root square and stopped 
the first time that the walk hits a square on the boundary of the quadrangulation). Then 
one can define a discrete scaling exponent (analogous to the box counting exponent in (J3J), 
with 5 replaced by n" 1 ) as follows: 



Identifying X n with its image in a conformal map to, say, D, one might guess that the random 
pair (X n , A n ) has a scaling limit (X, A), where X is a random subset of D (in our example, 
it might be a Brownian motion) and A is some form of the Gaussian free field. 

If this is the case, then on a heuristic level, one would expect that the quantum scaling 
exponent of X is A = A#, since, in the notation of Corollary 11.71 if we write 5 = n -1 , we 
would expect that E[5X(//, S, X)] scales like E(n _1 |X n |). 



7 2 = K = 8/3. 
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In discrete quantum gravity models, it is often possible to compute Ad explicitly (and 
rigorously) using random matrix techniques or tree bijections; it is also often possible to 
compute 7 directly using discrete quantum gravity machinery and so heuristically obtain its 
value in the continuum limit. 

Assuming values for A^> and 7 — and assuming A = A^ — the KPZ formula gives the 
Euclidean scaling dimension of X. In many interesting examples, X is a random fractal 
(a Schramm-Loewner evolution, for example, or the outer boundary of a planar Brownian 
motion) whose Euclidean scaling dimension might not be immediately obvious otherwise. 

Finally, we mention that, in the standard realm of conformal field theory, there exists a 
precise relation between the central charge c < 1 of the statistical model coupled to quantum 
gravity and the value of Liouville parameter, 7 = (V25 — c — y/l — c) / [ KPZ88tfDav88b 
IDK891 ISei90t IGM93] , as well as the corresponding connection between SLE K and Liouville 
quantum gravity models with 7 = min{/€, 16/k}. 

Our result extends the validity of the KPZ relation outside that CFT framework to any 
value of Liouville parameter 7 < 2, with the Ansatz that the fractal set X and the GFF are 
sampled independently. A possible interpretation of the KPZ relation in that case would 
be that it describes the quantum geometry of the given fractal in the quenched random 
metric generated by random graphs, equilibrated with a conformally invariant system with 
a value of c or k corresponding to the chosen value of 7. For example, one could first 
choose a random graph weighted by the critical Ising model partition function; and then 
perform a loop erased random walk on that graph, ignoring Ising clusters. In this case, one 
would expect the Euclidean dimension of the path to be that of SLE 2 (which corresponds to 
loop erased random walk), while the value of 7 describing the metric would be \/3 (which 
corresponds to the critical Ising model), and one could use KPZ to predict the quantum 
scaling dimension. 

Similar ideas appeared in previous numerical work ABT~99| IJJ99] , but the data so far 
appear as inconclusive. 

8 Future work 

The second author currently has two papers in preparation which aim to give additional 
support to the conjectures in the previous section. This section contains a brief outline of 
these forthcoming results. A joint paper from a physics perspective is also in preparation 
about the relation between Liouville quantum gravity and SLE. 

8.1 Quantum wedges and conformal welding 

In this section, we assume that < 7 < 2 has been fixed. Recall from Proposition 12.11 that 
two pairs (D, h) and (D, h) have the property that the quantum measure on one is the image 
of the quantum measure on the other under the conformal map ip : D — > D if and only if 

h = hoil) + Q\og\il)'\. (32) 

We define a "metric" to be an equivalence class of pairs (D, h) under the following 
relation: (D, h) ~ (D, h) if there exists a conformal map ip for which (J32l) holds. This 
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is a metric in the sense that areas and lengths can be computed via the definitions of 
this paper. We define a "metric with k marked points" to be an equivalence class of the 
set of triples (D, h, (z\, z 2 , . . . , Zk)) under the following relation: (D, h, (zi, z 2 , . . . , Zk)) ~ 
(D, h, (zi, z 2 , . . . , Zk)) if there exists a conformal map ip for which (1321) holds and ip(zi) = Zi 
for each 1 < i < k. 

Roughly speaking, a quantum wedge is the metric obtained by taking h to be an 
instance of the free boundary Gaussian free field on an infinite wedge {z : < arg z < 9} for 
some 9. If we conformally map the half plane to the wedge via the map ipe(z) = z 9 ^ n , then 
we see that this is also the metric defined by taking D — EI and letting h be an instance of 
the free boundary Gaussian free field on EI plus the deterministic function Qlog|^| which 
up to additive constant is Q log \ z e ^~ l \ = Q(9/n — 1) log \z\ = —a log \z\ for a := Q(l — 9/n). 
Our precise definition of quantum wedge (below) will include any aeK less than Q (which 
is the limit as the wedge angle tends to zero; note that our definition will allow for wedge 
angles greater than 2ir). 

Now, we have to be a bit careful since h is only defined up to additive constant, and 
adding a constant to h (thereby "rescaling the metric by a constant") does not yield an 
equivalent metric. However, if we define ifj b (z) = bz, for some b > 0, then 

h = ho + Q hg \(ilj b )'\ =hoip b + Q\ og b 

does yield an equivalent metric. Thus, rather than defining h modulo additive constants, 
we choose h modulo transformations of this form. (To be more precise, note that under the 
map ip b (z) the function B t = h e -t transforms as B t — > B t -i og b + Qlogb. In this case, the 
graph (t, B t ) in the t and B t plane is translated along the span of the vector (—1, Q). Two 
B t functions can arise from the same metric — and may be considered "equivalent" — if one 
is obtained from the other by a translation of the plane preserving the span of (—1, Q). We 
can define a canonical representative of the equivalence class by translating along the line 
so that B t first hits this line at the origin. In a quantum wedge, the law of B t for positive 
t is simply the law of Brownian motion with the given drift and Bt for negative t is the 
same but conditioned not to hit the line before 0. The law of the difference between h and 
its expectation given Bt turns out to be well defined without any undetermined additive 
constants.) 

The result is a random metric with two marked points (0 and oo) whose law is invariant 
under constant rescalings (Mobius transformations fixing those two points). Each quantum 
wedge — when parameterized by EI as discussed above — has an infinite amount of quantum 
mass, almost surely, but only a finite amount corresponding to any particular bounded subset 
of EI. (In particular, the law of a quantum wedge is not symmetric under reversing the two 
marked points, since every neighborhood of its second point has infinite mass, and this is not 
true of the first point.) The weight of the quantum wedge is the number defined from a as 
follows: W := 7(7 + ^ — 0;). Taking a less than Q corresponds to taking > 7(7 + ^ — Q) = 
7 2 /2. This definition is motivated by the following: 

Theorem 8.1. Choose a quantum wedge W of positive weight W , represented by some 
(D, h, (zi, z 2 ))- Suppose W = W\ + W 2 for some Wi > 7 2 /2 and then independently choose 
an SLE PliP2iK; for pi = Wi — 2 and k = 7 2 , from Z\ to z 2 . Let Y denote the set of points 
on the curve and let D\ and D 2 denote left and right components of D\T. (The condition 



35 



Wi > k/2 corresponds to the condition that the pi are in the range for which the path almost 
surely does not intersect the boundary ILSWOSf .) 

Then the random metrics W\ = (D\, h, (z%, z 2 )) (with h restricted to D\) and W 2 = 
(D 2 , h, (zx, z 2 )) (with h restricted to D 2 ) are independent. Each Wi has the law of a quantum 
wedge with weight Wi. 

In fact, one can use Theorem 18. II to give a definition of quantum wedges of positive weight 
less than 7 2 /2; these wedges are not wedges topologically (since their left side hits their right 
side at a random fractal set of points) but they can be well defined, and the above theorem 
holds for these wedges as well. It also turns out that W is uniquely determined by the 
Wi and may be obtained by conformal welding the right side of W\ to the left side of W 2 , 
where each is parameterized by quantum length. (This is closely related to a conjecture due 
to Peter Jones that SLE can be defined via conformal welding and the boundary measure 
induced by the Gaussian free field.) If we use the boundary analog of the correspondence 
in (|26|) to relate A and W, we find that this fact is equivalent to the additivity of quantum 
exponents predicted and advocated by the first author, which is well motivated in discrete 
quantum gravity models |Dup98 , Dup99b Dup04| . One can also weld the two sides of a 



wedge of weight W to each other to obtain a quantum cone (which is defined identically 
to the quantum wedge but with the whole plane instead of the half plane) with a' — f + -. 
We refer to this object as the quantum cone of weight W. 



8.2 Scaling limits of FK clusters on random graphs 

Consider a domain D with marked points Z\ and z 2 , and sample a pair (h, z) from (so 
that z is an interior marked point) and consider an instance h of the GFF on D with 
Dirichlet boundary conditions. Then choose a space-filling SLE K / path V from Z\ to z 2 , with 
k' = 16/ k = I6/7 2 . (When 4 < k' < 8, the path should be the space-filling analog of SLE K / 
constructed via the exploration trees corresponding to the conformal loop ensemble CLE K i, 
as advocated in [She].) 

Next, take an infinite volume limit by zooming in near the point z. The limiting object 
is a quantum cone with a = 7, together with a random space filling path that comes from 
00, hits z, and then fills up the rest of space. We parameterize the path r by M in such a 
way that T(0) = and T([a, b]) for each a < b in IR, has quantum measure b — a. 

It turns out that T[0, 00) (the "future") and r((— 00, 0]) (the "past") are independent 
quantum wedges. In fact, since the whole process is stationary with respect to t, this 
will imply an independence of increments result: namely that T([a, b}) and r([c, d]) are 
independent as random metrics when the (a, b) and (c, d) do not overlap. 

Denote by L t the quantum length of the left boundary of r((— 00, t}) minus the length 
of the left boundary of T((— 00, 0]) (both boundaries have infinite length; but since the two 
boundaries agree outside of some finite region, this can be defined as the difference in lengths 
of the portions that do not agree). Define Rt similarly. It turns out that the processes L t + Rt 
and L t — Rt are independent Brownian motions in the plane, but run at different speeds; for 
4 < k', we have 



Var[L, - Rt] 
Var[L t + Rt] 
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Note that this ratio ranges from oo to as k = 16/k' ranges from to 4. It is less than 1 
for k, G (2, 4) and equal to 1 for k = 2. 

We will also define an exploration path corresponding to a critical FK cluster model 
with parameter q. Given this relation, the analogously defined (L n , R n ) for these processes 
converge to the (L t , Rt) above under the usual Brownian rescalings. The variance ratios agree 
if we set q = 2 + 2cos^y = 2 + 2cos4r. This relationship between q and k' is consistent 
with a prediction made by the first author |Dup03j. In the case of k — 2, this convergence is 
an immediate consequence of a beautiful bijection by Bernardi, which shows that (L n , R n ) 
is a simple random walk on Z 2 [Ber06j. When q > 4, the metrics converge to the continuum 
random tree (a.k.a. branched polymer). 

The results above will show (barring technical difficulties), that the discrete quantum 
gravity loop models converge to a quantum cone endowed with a conformal loop ensemble, 
independently of the metric — however, this convergence is only in a very special topology, 
the topology in which metrics are considered close if the corresponding (L.,R.) processes 
(which we call driving functions) are close. Strengthening the topology to prove a re- 
sult like Conjecture 17.11 requires in some sense showing that when two processes have close 
driving functions, there is a high probability that they are close in the sense that, when 
both metrics are conformally embedded in the plane, the images of the exploration paths 
are close. This seems to reduce to proving a kind of fancy "random walk in random environ- 
ment" uniformization result, which, though intuitively plausible, appears extremely difficult 
to handle. Another way to strengthen the topology would be try to define a distance function 
on Liouville quantum gravity and prove convergence with respect to a topology on the set 
of metric spaces. This appears to be even more difficult, although the recent progress on 
topological scaling limits of discrete planar maps by Le Gall [LG07] and others may prove 
helpful in the special case k = 8/3, as well as the recent study of geodesies in large planar 
maps and in the Brownian map [LG08J (see also |BG08al BG08bj). 
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